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Abstract. We study the entropy production of a microscopic model for 
nonequilibrium wetting. We show that, in contrast to the equilibrium bound 
interface in a nonequilibrium steady state produces entropy. Interestingly, in some 
regions of the phase diagram a bound interface produces more entropy than a free 
interface. Moreover, by solving exactly a four-site system, we find that the first 
derivative of the entropy production with respect to the control parameter displays 
a discontinuity at the critical point of the wetting transition. 



1. Introduction 

For systems in thermal equilibrium the probability distribution of states is given by the 
Boltzmann-Gibbs measure, allowing one to calculate various macroscopic observables. 
Although such a theoretical concept is not available for systems out of equilibrium, it is 
nevertheless possible to make certain general statements in the nonequilibrium case. The 
most important statements of this kind are fluctuation relations [THT-4"] , which constrain 
the probability distribution of fluctuating entropy along a given stochastic trajectory of 
microstates. A less general result, analogous to the second law of thermodynamics, is 
that the average entropy production is non- negative [T5] . 

The average entropy production in the stationary state is a signature of 
nonequilibrium because it is zero if and only if detailed balance holds. Recent 
studies on the average entropy production explore its relation with the nonequilibrium 
steady state measure [T&H20] and its behavior at the critical point in systems with 
nonequilibrium phase transitions [5TH2S]. Interestingly, in rather different models the 
entropy production was found to peak near the critical point. Moreover, the first 
derivative of the entropy production with respect to the control parameter was found 
to diverge at the critical point for the majority-vote model [22] and in a nonequilibrium 
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Figure 1. Transition rates for the RSOS model. 



Ising model in contact with two reservoirs [24"||25]. However, it remains unclear whether 
these are general properties of nonequilibrium models with continuous phase transitions. 

In this paper we address these questions by studying nonequilibrium wetting 
transitions [26H28] which are known to exhibit a rich critical behavior out of equilibrium. 
Such transitions occur in models of growing interfaces belonging to the Kardar-Parisi- 
Zhang [29] universality class in the presence of a hard wall. The wetting transition is 
controlled by a parameter such as the growth rate, depending on which the interface 
either detaches or stays bounded to the wall. 

As a representative of this class of models we investigate the so-called restricted 
solid-on-solid (RSOS) model with a hard wall at zero height [30], extending a previous 
work where where we studied the entropy production of the free interface without a 
wall [18]. The presence of a hard wall changes the scenario completely: a bound phase 
arises in the region of the phase diagram where the free interface velocity would be 
negative. We show that a bound interface can produce entropy if detailed balance is 
not fulfilled and, moreover, in some regions of the phase diagram a stationary bound 
interface has a higher entropy production than a moving free interface. 

We also calculate the entropy production exactly for a four-site system. 
Interestingly, concerning entropy production, this small system displays the same 
qualitative behavior as an infinite one. As our main result, we show that in the non- 
equilibrium case, where detailed balance is violated, the first derivative of the entropy 
production displays a discontinuity with respect to the control parameter at the critical 
point. 

The paper is organized as follows. In the next section we first recall the definition 
of the model and briefly explain the wetting transition. In Sec. 3 we define entropy 
production, discuss its interpretation in the RSOS model and show numerical simulation 
results. Our main findings come in Sec. 4, where we present the exact calculations for a 
four-site system. Finally, after concluding in Sec. 5, some intuitive ideas about entropy 
production are discussed in the appendix. 
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2. The model and its phase diagram 

The RSOS model studied in this paper is defined on a one-dimensional lattice with L 
sites and periodic boundary conditions. The interface configuration is characterized by 
a height profile {h} = (hi, h 2 , . . . , hi,), where hi e Z is the height at site i. The interface 
evolves random-sequentially by spontaneous deposition and evaporation of particles with 
certain rates which are shown in Fig. [TJ Moreover, the dynamics is constrained by the 
restriction \h i+ i — hi\ — 0, 1, which introduces an effective surface tension. The model 
is controlled by two parameters, namely, the growth rate q and the rate for evaporation 
from the middle of plateaus p, while the evaporation rate from edges is set to 1. As 
usual, the probability Phi,...,h L if) to find the system at time t in the configuration {h} 
evolves according to the master equation. 

In this paper we are primarily interested in the interface velocity and the entropy 
production. These quantities can be expressed in terms of the three-site probability 
distribution 

oo 

PfH-uhiM+iit) = ^2 Ph u ...,h L (t) , (1) 

hi,...,hi-2,hi+2,—,hL=0 

which is obtained by integrating out all other height variables. For periodic boundary 
conditions this distribution is translationally invariant and thus independent of i. From 
the dynamical rules in Fig[H one can easily see that the interface velocity can be written 
in the form 

V = q ^~]{Ph,h,h + Ph+l,h,h + Ph,h,h+1 + Ph+l,h,h+l) 
h 

— P 2^ Ph,h,h — '5^ J (Ph-l,h,h + Pfi,h,h-1 + Ph-l,h,h-l) , (2) 
h h 

where we suppressed the argument t for the sake of readability. 

Depending on the control parameters q and p the interface of an infinite system 
will propagate (after an initial transient) with a constant positive or negative velocity v. 
This defines two regions in the phase diagram (see Fig. [2]), which are separated by the 
black line where v — 0. For given p we denote corresponding value of q as q c = q c {p)- 

Let us now assume that the growth process takes place on top of an inert substrate 
by adding a hard wall at zero height. Such a wall can be introduced by forbidding 
evaporation at height h — 0, meaning that negative heights become inaccessible [30lf3Tj . 
Obviously, this modification does not change the behavior above the black line q > q c , 
where the interface detaches from the wall and grows linearly in time with positive 
velocity. However, for q < q c the situation changes completely: Instead of a negative 
velocity, the interface is now pinned to the wall and thus its velocity is zero. Therefore, 
the model with a wall exhibits a wetting transition at q c from a bound phase for q < q c 
to a moving phase for q > q c . 
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Figure 2. Phase diagram of the RSOS model obtained from numerical simulations. 
Along the black line the asymptotic interface velocity is zero. The red line indicates 
the minimum of the free interface entropy production as a function of q for a given p. 

A natural order parameter of the wetting transition is the density of contact points, 
where the interface touches the wall. In terms of this order parameter, the wetting 
transition displays a very rich critical behavior, which has been the main focus of past 
research in the context of nonequilibrium wetting [261 - 128] . In the present study, however, 
we are not interested in the contact point density, instead we will focus on the behavior 
of the entropy production. 



3. Entropy production 

3.1. Definition of Entropy production 

A continuous-time Markov process, as taking place in the RSOS the model, is defined 
by a space Q of configurations c £ Q and rates w c ^ c i for spontaneous transitions from 
configuration c to configuration d . Denoting by P(c,t) the probability of being in 
configuration c at time t, the average entropy production in the environment (medium) 
is defined by [TO] US) (see Appendix) 

S(t) = ^ P(c, t)ww In (3) 

c,c' 

Obviously, this expression requires all transitions to be reversible, i.e. 

w c ^ c , ^ <(=> w d ^ c ^ , (4) 

meaning that the parameter p in the RSOS model has to be positive. Of particular 
interest in this paper will be the average entropy production in the stationary state 
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denoted by S, which is given by the above formula with the stationary measure 
P(c) = lim^oo P(c, t). Clearly, if detailed balance holds, i.e. P(c)u> c _> c / = P(c')w c >^ c , 
the entropy production is zero. Because of the second law, S is expected to be positive 
in a nonequilibrium stationary state. 

3.2. Entropy production in the RSOS model with a wall 

Using the three-site probability distribution ([I]) and ([3]), the mean entropy production 
for the RSOS model per site s = S/L can be expressed as 

s = In q^2(P h ,h,h + Ph+i,h,h + Ph,h,h+i) 

h 

,h,h-l -l,h,h-l) 

h 

+ \n^Y j P h>Kh -\n^y j P hAh . (5) 

PT p h ' 

This expression is valid for the free interface, where the sums run over all integer heights, 
as well as for the bound case, where the summation is restricted to h > 0. As we divided 
by the lattice size, s is an intensive quantity in L so that we can conveniently carry out 
the limit L — > oo. 

For p = 1 and q < 1 the stationary bound state of the RSOS model with a wall is 
given by a Boltzmann-Gibbs-like exponential measure [5TJ|[3Tj 

L 

P(h 1 ,h 2 ,...,h L ) cx X {h}]Jq hi , (6) 

i=i 

where x{h} is 1 if the heights are non-negative and the configuration obeys the RSOS 
constraint and zero otherwise. In this situation the dynamics is known to obey detailed 
balance and hence the entropy production is zero. For p = 1 and q > 1, however, the 
model is in the moving phase so that detailed balance is not fulfilled. In this case one 
can show that the entropy production is positive and linearly related to the velocity by 
s = vlnq. 

As shown in the Appendix, for p / 1, even in the bound phase, the system is out 
of equilibrium because different microscopic configurations are associated with different 
chemical potentials. Fig. [3] shows the entropy production of the RSOS model with and 
without the wall as a function of q for various values of p. As expected, both entropies 
become identical at the critical point and coincide in the moving phase. Keeping p 
fixed and varying q the entropy production assumes its minimum at a certain value 
q{s m i n ). As can be seen in Fig. [3] and Fig. [21 for p ^ 1 the location of this minimum 
is smaller than the critical point q c . Interestingly, for p = 4 the entropy production of 
the bound interface reaches a local maximum below the critical point. Moreover, there 
is an extended region where it is larger than the entropy production of a free interface. 
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Figure 3. Entropy production s of the RSOS model for p = 0.1 (left), p = 1 (middle), 
and p = 4 (right) with and without wall. The vertical dotted line indicates the 
transition point q c . The simulations were done with system size L = 1024. 



This means that the presence of the wall, which suppresses the flux of particles between 
the system and the reservoir, increases the entropy production in the environment. 

4. Exact results for small system size 

The left and the central panels of Fig. |3] suggest that the entropy production of the 
interface with the wall (black curves) varies smoothly at the critical point for p < 1 
while the right panel indicates a discontinuous derivative for p > 1. However, the 
numerical simulations are not accurate enough to confirm or disprove this conjecture. 

As a way out, let us consider finite systems. In contrast to usual phase transitions, 
which exist only in infinite systems, the RSOS model with a wall exhibits a wetting 
transition even if the system size is finite. Although the critical behavior of such a 
finite RSOS model is completely different from the usual critical behavior in the limit 
L — > oo j32j|33] , we observe that entropy production is remarkably robust and that even 
very small systems display essentially the same qualitative behavior of this quantity as 
infinite ones. 

The smallest non-trivial system size L = 3 is not appropriate for the present study 
because in this case detailed balance always hold in the bound phase. Therefore, we 
will focus on the case L — 4, which was solved exactly in Ref. [52] • In this case we have 
(apart from vertical shifts) 19 possible interface configurations. Because of translational 
invariance, only six of the 19 associated probabilities are independent, denoted as 

%h = Ph,h,h,h 

Uh — Ph,h,h,h+1 — Ph,h,h+l,h = Pfi,h+l,h,h = Ph+l,h,h,h 

z h-l — Ph-l,h,h,h — Ph,h-l,h,h = Pfi,h,h-l,h = Ph,h,h,h-1 

Uh = Ph,h+l,h+l,h = Ph,h,h+l,h+l = Ph+l,h,h,h+l = Ph+l,h+l,h,h 

v h — Ph,h+l,h,h+l = Ph+l,h,h+l,h 
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Figure 4. Exact phase diagram for L — 4. The dashed lines show the corresponding 
numerical data for L — > oo taken from Fig. [2J 



— Ph,h+l,h,h-l — Ph-l,h,h+l,h ~ Ph,h-l,h,h+l ~ Ph+l,h,h-l,h ■ (7) 

For a free interface the RSOS model can be mapped onto charges jumping on a ring 
by introducing charge variables <7; = h i+ i — hi = 0, ±1 [M]. Therefore, the problem 
of calculating the stationary measure of the RSOS model with L = 4 sites reduces to 
the problem of finding the eigenvector associated with the zero eigenvalue of a 6 x 6 
matrix [IB] . Following this procedure we can compute the entropy production §5§ and 
the velocity fl2]) for the free interface, 

s (free) / j = 16g 5 + 8q\p + 1) + 2q\q + p) (p - 1) - 2p(3 +p)(2q + l+ p) 

N(q,p) 

8q 5 + Aq\p + 1) - 2g 2 (g + p)(p - 1) - 2p(l + p)(2g + 1 + p) , g 
H In — 

iV(g,p) p 

r fr ee) f fl x _ 4(1 + g + 2g) (3g 4 - p(2 + p)) 

W ' PJ " iV(g,p) ' lJJ 

where the denominator is given by 

N(q, p) = 6p 3 + 10p 2 g 2 + 26p 2 g + 22p 2 + 48pg 3 + 59pg 2 (10) 



58pg + 17p + 47g 4 + 47g 3 + 27g 2 + llq + 2 



Setting ?/ frec ) (g, p) = gives the critical line 



p=-l + Vl + 3g c 4 . (11) 

In Fig. S]we plot the critical line together with the line where the entropy production 
([8]) is minimal. As in the case L — > oo, the line of minimal entropy production lies 
always slightly below the critical line and touches it at the equilibrium point p = 1. 

Let us now turn to the entropy production of a bound interface. As the wall 
breaks translational invariance in height direction, it is no longer possible to map the 
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Figure 5. Transition rates for L = 4 in terms of the six variables defined in Eq. (J7J. 



problem onto charges jumping on a ring. Therefore, we use a different method which 
was introduced in Ref. [32] to compute the density of contact points. Applied to the 
entropy production this method works as follows: In terms of the six variables ((Zj) the 
dynamic rules in a system with L = 4 sites form a network of transitions which is shown 
in Fig. |5j With the help of this figure it is easy to see that the corresponding master 
equation is given by 

x h = Ay h + Aqz h ~i -A[q+p(l- 5 h ,o)}x h 

Vh = q%h + 2it ft + v h + qw h - [1 + p(l - 6 h>0 ) + 3q]y h 

Zh = pxh+i + qv h + 2qu h + w h+1 - (p + 2 + 2q)z h 

il h = 2qy h + 2z h - 2(1 + q)u h 

v h = 2qy h + 2pz h - 2(p + 

u;^ = g^-i + py h - (q + l)w h , (12) 

where Wq = and all variables for negative heights are zero. Note that the six variables 
represent 19 different configurations, which is taken into account by certain multiplicities 
in the equations. For example, the factor 4 in term 4yh in the first equation comes from 
the fact that there are four different height profiles represented by y^ and they all go 
to Xh by means of an evaporation event at rate 1. Moreover, note that translational 
invariance in height direction is broken by the height-dependent terms on the right-hand 
side. For example, the terms (1 — 5h,o) come from the fact that in the presence of a wall 
evaporation is forbidden at zero height. 

The stationary solution of the master equation ( [T2l) reads 



7 -i ( 3g 4 \h 3q 3 
Xh = Z ~77T I — \ ' Vh = Q x h, z h = — — x h , w h = px h , 
\p{2 + p)/ 2 + p 

( q 2 3g 3 \ / q 2 3pq 3 \ 

Uh= t—l + (q+l)(2 + p)) Xh ' Vh= {—l + (q + m +P )) Xh ' (13) 
where Z is a normalization factor. This factor can be determined by computing the 
normalization YlhLo = 1 of the single-site probability distribution Ph, which in terms 
of the six variables is given by 

Ph = x h + 3y h + Uh-i + z h + 3z h ^i + 2u h 

+ 2u h - 1 + v h + v h -i + 2w h + w h -i + w h +i (14) 
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Figure 6. The entropy production ([TT]) s(q) for p = 0.1 (left) and p = 4 (right). 



SO that 



r-l 



2p + p 2 - 3q 4 



p(2 + p + 8q + Apq + 12q 2 + 6pq 2 + 12q 3 + 12q 4 



(15) 



Finally, expressing the three-site probability distribution ([T]) in terms of the four-site 
probability distribution ([7]) and using formula fl5]), we obtain 



n (bound) 



(?,P) 



2(p- i)<r 



1 + g)(2 + p + 8q + Apq + 12q 2 + 6pg 2 + 12q 3 + 12q 4 



hip . 



(16) 



which is valid in the stationary bound state below the critical line ( ITTj) . Therefore, the 
entropy production of the model with a wall at zero height can be expressed as 



s(q,p) = s^ d \q,p)9(q - q c ) + s^ cc \q,p)9(q c - q) 



(17) 
and 



where 9(x) is the Heaviside step function, q c is given by (HU), s (frcc) (g,p) by i 
s {hounA \q,p) by (ITBI) . 

In Fig. [6] we plot the entropy production (Tl7j) as a function of q for two different 
values of p. It is clear that it is continuous at the critical point, which can be verified 
analytically by checking that ([8]) and ffTBI) become identical along the critical line. 
Contrarily, the first derivative with respect to q, denoted as 



d 

C(q,p) = —s(q,p). 



;is) 



is found to be discontinuous at the critical point for p ^ 1. The explicit expression for 
this derivative is too cumbersome to be written down. However, plotting the derivative 
in Fig. [7] for p = 0.1 and p = 4 one can clearly see a jump at the critical point. This 
is different from the results obtained in other nonequilibrium models [2211241125] . where 
the first derivative was found to diverge at criticality. 

The only exception is the case of detailed balance p = 1. Here the entropy 
production of the free interface varies as \nq(q — 1) and vanishes at the critical point 
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Figure 7. The first derivative of the entropy production C(q,p) for p = 0.1 (left) and 
p = 4 (right). 



q c — 1 while the entropy production of a bound interface vanishes for any q < q c . 
Therefore, the first derivative is continuous for p = 1, independent of the system size. 
We therefore conclude that the observed discontinuity is a characteristic property of 
nonequilibrium in the present model. Morevoer, we conjecture that the size of the jump 
is a measure of how far the system is driven away from equilibrium. 



5. Conclusion 

In this paper we have studied entropy production in the RSOS model with a wall which 
displays a nonequilibrium wetting transition. We have shown that in nonequilibrium 
a bound interface produces entropy due to an imbalance in the chemical potential 
associated with different local microscopic configurations. Surprisingly, we find that 
in some regions of the phase diagram a nonequilibrium bound interface produces more 
entropy than a moving one. 

We have solved the model exactly for the case L — 4, allowing for a precise analysis 
of the entropy production at the critical point. The main result of the paper is that 
the first derivative of the entropy production with respect to the control parameter q 
displays a discontinuity at the critical point. Furthermore, such a discontinuity occurs 
only under nonequilibrium conditions, whereas for p — 1, where detailed balance is 
fulfilled in the bound phase, the first derivative is found to be continuous. 

The entropy production was found to peak near the critical point in other 
nonequilibrium phase transitions [21T425] and its first derivative was found to diverge 
at the critical point [221 [2U • In this paper we showed that in nonequilibrium 
wetting transitions its first derivative displays a discontinuity at criticality. It would 
be interesting to further investigate the critical behavior of the entropy production in 
other systems out of equilibrium, it might be that the entropy production plays an 
important role in classifying different nonequilibrium phase transitions. 
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Appendix A. Motivation of the formula for entropy production 

Formula ([3]) gives a mathematical quantity that, in principle, can be calculated for any 
Markov process. The name entropy production of the external medium can be justified 
as follows. Consider a system (S) in contact with a thermal reservoir (R) at temperature 
T with the stationary measure P(c) oc exp[— E(c)/T]. For such a system the detailed 
balance condition reads 

^— = -hm-E(c))- (a.i) 

If the stationary system jumps from configuration c — > c', its energy changes by 
AE$(c — > c') = E(c')—E(c) while the energy of the external reservoir changes oppositely 
by AE R (c — > c') = — AE s (c — > c'). As the temperature of the reservoir is constant, the 
detailed balance condition (lA.ip implies that 

Ag= A ^ (c ^ c/) =ln^. (A.2) 

T w c >^ c 

In order to make the formula for entropy production plausible, let us now consider a 
fictitious system which is randomly switching between two thermal reservoirs at different 
temperatures T\ and T2. Such situation is described by a Markov process with transition 
rates w c >^ c = w^\ c + w^\ c , where the rates w^\ c and w^\ c are the equilibrium rates 
for each of the two reservoirs. 

In general, such a system is out of equilibrium and its stationary measure is 
not known. Let us assume that the system is in its nonequilibrium stationary state, 
undergoing a sequence of N transitions c —> c\ —>•...—>• cjv- Using formula (IA.2[) the 
variation of the energy in the reservoirs u = 1,2 can be written as 



N (ui) 

A^=T^^ln^p, 

i=l Wa^Ci-x 



where ^ = 1,2 labels the reservoir that the system was connected to during the jump 
Ci-i — > Ci. The total change of the entropy in the external medium (consisting of the two 
reservoirs) is then given by the sum of the energy differences in each reservoir divided 
by its respective temperature, i.e. 

A Ul) AF (2) N l Vi ) 

AS= ^ + ^ = Y>te. (A.4) 



Ci-1 



Finally, fixing a time interval At and averaging AS/ At over of all possible trajectories 
we obtain the formula for entropy production Q in the stationary state. Note that 
these arguments rely on the assumption that the two reservoirs are almost equilibrated, 
which is reflected by using (I A . 2 [) in (IA.3[) . This assumption is discussed in more detail 
in Ref. 135 1. 
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The above argument can be used to explain the positive entropy production of 
the RSOS model in the stationary state for p / 1 To this end let us consider the 
growth process as a system in contact with a reservoir of particles, where the chemical 
potential for the deposition and evaporation of particles depends on the local microscopic 
configuration. If the interface separates two different phases A and B, deposition and 
evaporation events can be interpreted as reactions B — >■ A and A — >■ B, respectively. 
Let us take the equilibrium case p = 1 and denote by fi^ an d \Ib the chemical potentials 
of A and B, respectively. Hence if N particles are deposited, the energy in the external 
reservoir will change by AEr = NAfi = N(hb — Ha), where q = e AM . 

Now let us assume that the system is coupled to two reservoirs, namely one reservoir 
with chemical potential difference A/ii for deposition and evaporation on flat parts 
and at the edges and a second reservoir for evaporation from plateaus with chemical 
potential difference A/i 2 . Assuming the reservoirs to be almost equilibrated we may use 
the detailed balance condition (lA.lj) to identify q = e Af11 and q/p = e AM2 . When the two 
chemical potentials are different (p ^ 1) the system is out of equilibrium and there is an 
average flux of energy between the system and the reservoir even in the bound phase. 
The flux components (growth velocities) related to the chemical potential differences 
A/ii and A/i 2 are 

Vl = q y~](-Pfe,ft,fe + Ph+l,h,h + Ph,h,h+l — Ph-l,h,h ~ Ph,h,h-1 ~ Ph-l,h,h-l)i 
h 

V2 = 1 P h+l,h,h+l -P/; Ph,h,h- (A. 5) 

h h 

Therefore, v = V\ + v 2 and s = A/x^i + A/i 2 ^ 2 . In the bound phase we have V\ = — v% 
and the entropy production is s = vi(Ap,i — AyU 2 ). Hence, out of equilibrium, the 
imbalance of the chemical potentials for different microscopic configurations causes the 
entropy production in the bound phase to be positive. 
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